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We investigate the decays of Xc0,2 into vector meson pairs (VV), pseudoscalar pairs (PP), and 
scalar pairs (5'5') in a general factorization scheme. The purpose is to clarify the role played by the 
OZI-rule violations and SU(3) flavour breakings in the decay transitions, and their correlations with 
the final state meson wavefunctions. For Xc0,2 —> VV and PP, we obtain an overall self-contained 
description of the experimental data. Applying this factorization to XcO —> /o/o, where i, j = 1, 2, 3 
denotes /o(1710), /o(1500) and /o(1370), respectively, we find that specific patterns will arise from 
the model predictions for the decay branching ratio magnitudes, and useful information about the 
structure of those three scalars can be abstracted. 
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Charmonium decays provide a great opportunity for studying the gluon dynamics at relatively low energies. In line 
with the lattice QCD predictions for the existence of the lowest glueball state (0'^+) at 1.45~ 1.75 GeV (i, one 
of the most important motivations to study charmonium decays is to search for glueball states and fill in the missing 
link of QCD. During the past few years, enriched information from experiments reveals a rather crowded scalar meson 
spectrum at 1^ 2 GeV. For example, three /o states, /q'^''^ (= /o(1710), /o(1500), /o(1370)}, are observed in both 
^sj PP scattering and pp annihilations with different decay modes into pseudoscalar pairs ^ Q, 0, 0|- They are also 
^ confirmed by the BES data in J/^|J ujfl and 0/g 0,IEI3- Moreover, another scalar /o(1790), which is distinguished 
Q . from /o(1710), is also seen at BES 0. In contrast with the /o(1710), /o(1790) does not obviously couple to KK, 
' while its coupling to tttt is quite strong. Since all these states cannot be simply accommodated into a QQ meson 
^ , configuration, observation of such a rich scalar-meson spectrum could be signals for exotic states beyond simple QQ 
' configurations The fast-growing database also allows us to develop QCD phcnomcnologies from which we 

can gain insights into the glueball production mechanism 0, Q, ITEL ITfiL Il7l | . 
fSj ' So far, it is still a mystery that the glueball states have remained hidden from experiments for such a long time 
J> ' 01 . One possibility could be that a pure glueball state cannot survive a long enough lifetime before it transits into 
^ ■ QQ pairs. Because of this, a glueball would be more likely to mix with nearby QQ. As a result, it should be more 
' sensible to look for a glueball-QQ mixing state, instead of a pure glueball. Such a mixing scheme, on the one hand, 
■ indeed highlights some characters of the meson properties, and fits in the experimental data quite well 0,0; 01- It 
' seems also promising in accounting for the crowded scalar meson spectrum at 1~ 2 GeV [l^. On the other hand, 
ly-j , it raises questions about the role played by the OZI rule |0 in the glueball production mechanism, which needs 
' understanding in a much broader context. In Ref. "231, it is shown that the idea of intermediate virtual meson 
exchanges (2ll |2^ |2^ can provide a dynamic explanation for the large OZI-rule violations in J/i/; Li;/o(1710) and 
(^/o(1710). It implies that correlations of OZI-rule violations could be an important dynamic process in charmonium 
I ' decays. Such correlations should not be exclusively restricted to the scalar meson productions in J /ip decays. It 
could also contribute to other processes ,2^], such as Xc0,2 hadronic decays. Therefore, a systematic investigation of 
^ XcO, 2 — > VV, PP and SS is necessarily useful for clarifying the role played by the OZI-rule violations. 
. • , The other relevant issue in the charmonium decays into light hadrons is the SU(3) flavour symmetry breakings. In 
many cases, the flavour-blind assumption turns out to be a good approximation for the gluon-QQ couplings. However, 
due to the presence of OZI-violations, an explicit consideration of the SU(3) flavour symmetry breaking will be useful 
for disentangling correlations from different mechanisms. 
d In this work, we will present a factorization scheme for Xc0,2 — > VV , PP, and SS. Our purpose is to disentangle the 

roles played the OZI-rule violations and SU(3) flavour symmetry breakings, which will correlate with the final state 
meson wavefunctions in Xc0,2 decays. By this study, we hope to abstract additional information about the structure 
of the scalar mesons /q. It can be reflected by specific patterns arising from the decay branching ratio magnitudes 
due to such correlations. 

To proceed, we first distinguish the singly OZI disconnected (SOZI) processes and doubly OZI disconnected (DOZI) 
ones by considering transitions illustrated by Fig. ^ The SOZI transition [Fig. H^a)] hence can be expressed as 

(gi92g3g4|Vb|g5) = 5(i3>5{24) , (l) 

where Vb is the interaction potential, and \gg) denotes the minimum two-gluon radiations in cc annihilation; 1713 and 
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(724 are the coupling strengths. For the non-strange couphng, we denote 313 = (724 — 9o- To include the SU(3) flavour 
symmetry breaking effects, we introduce R = {ss\Vo\g) / {uu\Vo\g) = {ss\Vo\g) / {dd\Vo\g) , for which i? = 1 refers to the 
SU(3) flavour symmetry limit. 

The DOZI process of Fig. ^h) describes the transition that the two qq pairs are created and recoil without quark 
exchanges. The amplitude can be expressed as 

{qiq2q3q4\Vi\gg) = 5(i2)5{34) = r9(i2).g(34> , (2) 

where Vi denotes the interaction potential, and parameter r is defined as the relative strength between the SOZI and 
DOZI transition amplitudes. It is introduced to take into account the DOZI effects. 

The two gluons from the cc annihilation can also couple to the glueball components via Fig. . The transition 
amplitude can be expressed as 

{GG\V2\gg)=g"g" = tg^o . (3) 

where V2 is the interaction potential and g" is the gluon coupling to the scalar glueball; Parameter t is introduced to 
account for its coupling strength relative to the SOZI process. 

Another possible process in the decay of Xc0,2 is via Fig. ^d), where the two gluons will couple to qq and glueball, 
respectively. According to Eqs.ElandEl its relative coupling strength to the SOZI process can be written as 

{qiq2G\V3\gg) = V^g^ . (4) 

With the above factorization scheme, we consider a transition of Xc0,2 M'^M^ , where and denote two 
final-state quarkonia such as VV, PP, or /q/q- Firstly, in a general basis we assume a wavefunction for / = M*: 

AP = x^\G) + y,\ss) + z,\nn) , (5) 

where \G), |ss), and \nn) — \uu + dd)/\/2 are the pure glueball, and flavor singlet qq components, respectively; 
Coefficients Xi, yt and Zi denote the mixing angles of those components, and satisfy the unitarity condition. 
The transition amplitude for XcQ,2 M'^M^ can be written as 

{APA-P\{Vo + Vi + V2 + V^)\gg) 
= {{xiG + yiss + z.inn)(xjG + Vjss + z.jnn)\{yo + Vi + V2 + V:i)\gg) . (6) 

Notice that the potential Vo will only allow the transitions to / = meson pairs, i.e. transitions of {nnss\Vi\gg) — 
{ssnn\Vi\gg) = 0. We can then decompose the above transition into 

{KPAP\{Vo + Vi + V2 + V3)\gg) 
= gQ[xi{tXj -t- RVrtyj + \j2rtzj) -|- Ryi{\/rixj + (1 + r)R.yj + \pirzj) 

+Zi{V2Hxj + V2rRy.j + {1 + 2r)zj)] . (7) 

For Xc0,2 decays into other octet meson pairs with / 0, e.g. KK, tttt, etc, we can see that only Vq transition is 
allowed, and all the others are forbidden. This selection rule immediately requires that decay channels with 1^0 
will not suffer from the OZI-rule violations. As a consequence, the OZI violation effects and SU(3) flavour symmetry 
breaking can be separated out. The explicit reduction of Eq. Qfor XcO.2 VV and PP will be given later. 

In the calculation of the partial decay widths, a commonly useful form factor is applied: 

T^{p)^p''exp{-py8(3') , (8) 

where p and I are the three momentum and relative angular momentum of the final-state mesons, respectively, in the 
XcO.2 rest frame. We adopt f3 = 0.5 GeV, which is the same as in Refs. [Ulllll^. Such a form factor will largely 
account for the size effects from the spatial wavefunctions of the initial and final state mesons. After taking into 
account this, we assume that parameter 170 will be an overall factor for a given initial charmonium. 

i) Xco,2 ^VV 

For Xc0,2 (t>4': ^<^7 a-iid LU(j), the transition amplitude reduces to simple forms due to their ideally mixed flavor 
wavefunctions, i.e. (j) = ss, and w = nn = {uu + dd)/y/2. This means that for (j) meson, Xi — Zi = and yi = 1, and 
for Lo meson, Xj — yj — and Zi = I. Therefore, we have 

mV\99) = 5o^'(l + 
{uu;\V\gg) = 5o(l + 2r) 

(iomgg) - g',rRV2 , (9) 



where we compactly write Vq + Vi +V2 + V3 as V". 
For Xc0,2 — K* K* and pp, we have 



{K*+K*-\V\gg) = g^R^ 

{p^p-\V\gg) = gl ■ (10) 

It can be shown that the amphtudes for other charge combinations are the same. 

Interesting correlations arise from those decay amplitudes. Equation 1101 shows that the decay of Xc0,2 PP is free 
of interferences from the DOZI processes. Nevertheless, it does not suffer from the possible SU(3) flavour symmetry 
breaking. Ideally, a measurement of this branching ratio will be useful for us to determine go. Unfortunately, 
experimental data for this channel are not available. The decay of Xc0,2 K*K* also merits great advantages. It is 
also free of DOZI interferences, and the SU(3) breaking effects have simple correlation with the basic amplitude g^. 
In contrast, Eq. El involves correlations from both OZI-rule violations and SU(3) breakings. 

So far, three channels have been measured for Xco decays by BES collaboration |2Ml2al23: -B-^xco-"t>4> — (1-0 ± 
0.6) X 10-3, BR^^„^^^ = (2.29 ±0.58 ±0.41) x lO^^, and BR^^^^^,ok'0 = (1.78 ± 0.34 ± 0.34) x lo'^^^. This allows 
us to determine the parameters with very small x^- 

r = 0.45 ±0.48, i? = 0.90 ± 0.22, go = 0.25 ± 0.06 GeV^/^ (11) 

We note that the central values can be obtained by directly solving the equations for these three channels. Fitting 
the data can provide an estimate of the parameter uncertainties due to the experimental error bars. The results show 
that large errors of the data can lead to large uncertainties for the parameters, especially for r. To see more clearly 
their correlations, we take the fraction between the branching ratios for ujui and (fxj) to derive: 

_ R^Cq - 1 

2 - i?2Co ' ^^^> 

where Co = [p^jBR^^^^^^T"^ [Tpcj,) /PuiBR^^g^^^T'^ {tP^)]^/"^ = 1.6. It shows that for a range of 1/Co < R^ < 2/Co, 
the above relation leads to < r < +00. In other words, the combined uncertainties from the SU(3) breaking and 
experimental uncertainties can result in rather significant changes to r. Therefore, improved measurements are still 
needed for the determination of r. Meanwhile, the value for R suggests that the SU(3) breaking is about 20%. This 
seems to be consistent with the success of SU(3) flavour symmetry in many circumstances. 

Adopting the central values for the parameters, we can predict the branching ratios for Xco ~^ ^0 and pp, and the 
results are listed in Tabled We find a relatively smaller branching ratio for Xco — > ^4>^ i-e. BR^^^^^^ = 0.45 x 10^^. 
This channel should be forbidden if the OZI-rule is respected. However, with the errors from Eq. the estimate of 
the root mean square error gives 1.5 x 10"'^, which is the same order of magnitude as (pcf) and louj channels. This implies 
that the OZI-rule violations could be significant. A direct measurement of ujcj) channel will put a better constraint on 
the OZI-rule violation effects. 

The decay of Xco PP does not suffer the OZI rule, and its branching ratio is predicted to be sizeable. Notice that 
Xco has large branching ratios for 2(7r+7r~) and p^tt+tt" |2^, the predicted value, BR^^^^pp = 1.88 x 10"^, appears 
to be reasonable. A precise measurement of this channel is strongly desired for the determination of parameter gg. 
Note that go appears in the branching ratio with a power of four. The error for go as in Eg. Illl can still lead to large 
uncertainties of 92% for BR^^„-,pp. 

For Xc2 — ^ yy , so far, branching ratios for three channels have been available from BES [25ll26Ll2^ : BR^^^^^^ — 
(2.00±0.55±0.61) X 10-3, BR^^^^^^ = (1.77± 0.47± 0.36) x 10"^ and BR^^^_j^,j^, = (4.86 ± 0.56 ± 0.88)''x IQ-^. 
This will allow us to make a parallel analysis as the XcO decays. By fitting the data, we obtain 

r = 0.24 ±0.29, i? = 1.09 ± 0.21, go = 0.26 ± 0.06 GeV^/^ (13) 

Again, the central values can be obtained by solving the corresponding equations, and large uncertainties for the 
fitting results highlight the effects from the experimental errors. It is interesting that the parameters for XcO and Xc2 
are not dramatically different from each other. Instead, they are quite consistent. In particular, in contrast with XcO 
decays, the calculation shows that the OZI-rule is better respected in Xc2 decays. We also note that with a factor 3 
reduction of the data errors in both cases, the uncertainty of r can be improved to the second decimal place. 

With the central values of the fitted parameters, we predict the branching ratios for wt/) and pp channels. Similar 
to XcO: the channel tufj) is relatively suppressed with BR^^^—n^^tp = 0.24 x 10-^. The branching ratio for pp is found 
sizeable with a value of BR^^^-^pp — 2.41 x 10-^. We also present the estimated errors for these two channels in 
Table D 

ii) Xco,2 ^ PP 



The same analysis can be applied to the XcO.2 decays into pseudoscalar meson pairs. We will adopt the same 
notations as Part (i) for the parameters. But one should keep in mind that they do not necessarily have the same 
values. For / = channels, i.e. Xc0,2 VVi v'v' ^ rjrj' , the mixing of ij and rj' are defined as 



77 = cos ann — sm ass 

r] = sin ann + cos ass , (14) 

where a = 54.7° +9p and 6p is the octet-singlet mixing angle in the SU(3) flavour basis. We do not consider possible 
glueball components mixed within the rj and 77' wavefunctions at this moment though they can be included. Quite 
directly, the transition amplitude of Eq. Qcan be reduced by requiring Xi — Xj = and t = 0: 

{r]r]\V\gg) = 5q[cos^ a + i?^ sin^ a + r(V2cosa - i?sina)^] 

{ri'ri'lVlgg) = ^^[sin^ a + i?^ cos^ a + r(\/2 sina + i? cos a)^] 

{m'\V\99) = glr[{l- R^/2)sm2a + V2Rcos2a] 

{K+K-\V\gg) = g^R 

{7r+7T-\V\gg) = gl , (15) 

where the transition potential is a compact form for Vq + V\ + V2 + V3. 

Experimental data for 7777, K^K~ and rnr are available |2M . l2it l3fll | . which will allow us to constrain the parameters 
and examine the model. We first determine in Xc0,2 tttt. Then, by taking the ratio between K^K~ and tttt, we 
determine the SU(3) breaking parameter R: 



R 



TTTT 



(16) 



where factors 2 and 3 are the weighting factors for the charged kaon and pion decay channels. Substituting r^+^- = 
6.0 X 10~'^rtot and Ftttt = 7.4 x W~^Ttot H^HS into the above equation, we have R — 1.06, which is in good agreement 
with the range found for Xco ~* VV. 

With the data for Xco ~^ we, in principle, can determine r via the first equation in Eg. 1151 However, we find that 
within the large uncertainties of the present data, BR^^^^rm = (2.1 ± 1.1) x 10~^, the determination of parameter r 
also possesses large uncertainties. To show this problem, we plot the branching ratios for xm ~> v'v' ^-nd 777;' in 
terms of r in Fig. Efa). The solid line is for BR^^g^rjrj, which exhibits a slow change with r. The arrow acrosses the 
central value of the experimental data, and marks the error bars. Therefore, the slow variation of B R-^^^^rji) makes 
it difficult to determine r. Interestingly, it shows that BR^^^^rii^i' and BR^^^^riri' are both fast changing functions 
of r. In other words, one can determine parameter r by measuring these two branching ratios, and matching the 
pattern of the curves in a narrow region of r. Meanwhile, an improved measurement of 7/7; channel is also strongly 
recommended. 

For Xc2 PP, there are also experimental data available for three channels: BR^^^^k+ j^- — (9.4 ± 2.1) x 10"'*, 
-Bi?x^2_»7r+7r- — (1.77 ± 0.27) X 10~^, and BR^^^^vv < ^-^ ^ 10"'^. However, notice that the branching ratio for 7777 
is only an upper limit and quite large uncertainties are with the data, the parameters cannot be well-constrained. 
Similar to case of XcO PP, we can determine R and 50 with the data for tttt and KK: R = 0.70, go = 0.30. In 
comparison with the previous cases, it shows large SU(3) symmetry breakings. A possible reason could be due to the 
poor status of the data. 

Similar to XcO decays, we plot the change of branching ratios of 7777, rj'rj' in terms of a range of r. The results are 
presented in Fig.|21[b). Taking the upper limit of 7777 branching ratio as a constraint, the OZI-rule violation parameter 
r at least has a value of —0.16. Note that the relative branching ratio magnitudes among these three channels vary 
fast. Additional experimental information about 7777' or rj'rj' will provide better constraint on r. 

iii) XcO, 2 SS 

For XcO, 2 SS, where S denotes scalars, a direct gg GG coupling will also contribute to the 7 = decay 
amplitudes, and an additional parameter t has to be included (see Eq. ITJ. Therefore, new features may arise from 
the 1 = decay channels. In particular, it should be interesting to investigate the sensitivity of the decay branching 
ratios to the structure of the 1 = scalars, i.e. /q. Specific patterns arising from these channels will be useful for 
providing evidence for the presence of glueball components in the scalar meson wavefunctions. 

In contrast, the decays of Xc0,2 ^ ao(1450)ao(1450), and i^o(1430)i?o (1430) will be free of the interferences from 
the gg — *■ GG transitions due to the isospin conservation. 

To proceed, we will focus on the decays of XcO ^ /o/o- principle, all the decay channels into any two of 
these scalars are allowed except for /q/q, of which the threshold is slightly above the XcO mass. We will adopt the 



wavefunctions based on the glueball-QQ mixings [Ta . Il5l | , which are examined in the J/^ hadronic decays ^3 • As 
shown in Ref. T^, the branching ratios for J/^ 7f JTxhibit specific patterns, which cannot be explained by"mple 
QQ mixings within the scalars. It is also found in Ref. Jjgj that strong OZI-rule violation effects could contribute to 
the transition amplitudes. 

Due to lack of data, the parameters cannot be explicitly determined. However, based on our experiences from the 
first two parts, we can assume that the SU(3) flavour symmetry breaking is negligible. Thus, we fix i? = 1 in the 
calculations as an approximation. Also, for gg — > GG coupling, we assume that it has the same strength as the SOZI 
transitions, i.e. t = \. We are then left with parameters go and r for which, in principle, experimental constraints 
are needed. But we can still proceed one step further to consider the branching ratio fractions, where parameter g^ 
cancels. 

Taking Xco — > /o/o reference, we plot in Fig.|3|the branching ratio fractions of BRp ^ / BRfij?, in terms of r in 
a range of < r < 3. Interestingly, it shows that the branching ratio fractions are separated into two distinguished 
regions. In the limit of r ^ 0, the branching ratio of Xco ~^ fofo becomes the smallest one, and Xco ~* /o/o the 
largest. Howe ver at the region of r > 1, the branching ratio of XcO /o/o becomes the largest. 

In Refs. [la) 113' large OZI-rule violations are found in association with the scalar meson production in J/i/; decays. 
Therefore, we conjecture that large contributions from the DOZI processes will be present in XcO /o/o • Based on 
this, the branching ratio fractions at r > 1 strongly suggest that XcO /o/o '^i^l t)6 the largest decay channel. Notice 
that /o couples strongly to KK, while /q strongly to tttt, one might be able to see rather clear evidence for these two 
states in e.g. XcO ^ /o/o ^ K+K-n+n- at BBS [p. 

Due to lack of data for Xc2 — *■ SS, we do not discuss this channel in this work. 

To summarize, by distinguishing the SOZI and DOZI processes and including a direct glueball production mechanism 
in the transition amplitudes, the overall available data can be accounted for. Since these basic transition processes 
have been factorizcd out, we can study their correlations among all the decay channels. We find that the OZI-rule 
violations play quite different roles in the production of VV, PP, and SS in the Xc0,2 decays. This may explain 
that pQCD approaches systematically underestimate the data 1^3) 13 El] ■ On the other hand, a recent study 
based on the ^Pq quark pair creation model shows that the nonperturbative mechanism still plays important roles in 
XcO ~^ 4'4' US • In this approach, a coherent description of both pQCD favored and unfavored mechanisms allows us 
not only to isolate their contributions, but also to investigate their correlated interferences. Such correlations produce 
specific patterns for the branching ratio magnitudes of different channels, and can be highlighted by experimental 
data. 

For XcO, 2 ^ we find that the branching ratio of XcO ^ '-^4' is crucial for determine the role of the OZI-rule 
violations, and for Xco.2 PPi an improved measurement of will clarify the correlations from the OZI 

violations. For XcO ^ /o/o^ find that large OZI-rule violations will lead to suppressions on most of the branching 
ratios except /q/q = /o(1710)/o(1370). With an analogue to the large OZI-rule violations in the /o scalar productions 
'm J/tp ^ Vfo, we would expect that large OZI-rule violations would occur in XcO — > /o/o- Therefore, an observation 
of sizeable branching ratio for /q/q — /o(1710)/o(1370) will be a strong signal for the glueball-QQ mixings within 
the scalars. In contrast, small OZI violations will lead to relatively smaller branching ratios for /q/q . If this occurs, 
we will then have difficulty to understand the large OZI violations in J/-0 ^ V/q 0, |23| ■ Experimental data from 
BES Collaboration can provide an important test of this approach, and also put a more stringent constraint on the 
model parameters. We expect that a confirmation of this will also provide an additional evidence for the glueball-QQ 
mixings within those scalars. 

The author thanks F.E Close and B.S. Zou for useful comments on this work, and Z.J. Guo and C.Z. Yuan for useful 
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Energy, Chinese Academy of Sciences. 
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Decay channels 




Bi?xc2->vv(xl0"^) 


H 


1.0 ±0.6 


2.00 ±0.82 


U!OJ 


2.29 ±0.71 


1.77 ±0.59 




1.78 ±0.48 


4.86 ± 1.04 


ojcj) 


0.45 (1.07) 


0.24 (0.65) 


pp 


1.88 (1.80) 


2.41 (2.22) 



TABLE I: The branching ratios for Xc0,2 — > VV . The first three rows (00, uto and K* K*) are experimental data from 
BES (illilE^I, while the last two rows, highlighted by bold faces, are the model predictions for uo(j> and pp channels. The 
numbers in the brackets are the root mean square errors estimated with Egs. Illl and ll3l respectively. 




FIG. 1: Schematic pictures for the decays of Xc0,2 into meson pairs via the production of different components. 




FIG. 2: Branching ratios for the decays of XcO (Fig. (a)) and Xc2 (Fig.(b)) into r^rj (sohd curve), ri'rj' (dashed curve), and rjrj' 
(dotted curve) in terms of the SU(3) parameter r. The two-direction arrow in (a) denotes the central value of BR^^^^^ = (2.1± 
1.1) X 10~^, while the length marks the error bars. The thin straight line in (b) labels the upper limit BR^J^^^y^ < 1.5x10"^ jiffll ). 




FIG. 3: Branching ratio fractions for XcO decays into /o pairs in terms of the OZI-rule violation parameter r. The thin straight 
Una denotes ratio unity and separates the two regions for the ratios greater/smaller than one. 



